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Abstract 

We study optimal double helices with straight axes (or the fattest tubes around them) computa- 
tionally using three kinds of functionals; ideal ones using ropelength, best volume packing ones, and 
energy minimizers using two one-parameter families of interaction energies between two strands of 
types r~" and ^ exp( — fer). We compare the numerical results with experimental data of DNA. 
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1 Introduction 

In this papei0 we study three kinds of functionals for a one-parameter family of infinitely long double 
helices with straight axes that have radius 1 and slope a > 0, i.e. pitch P = 2T:a: 

r(a) =ri(a)ur2(a) (a>0), 
Ti{a) = {(cos 6*, sin 6*, 06*)}, r2(a) = {(- cos 6I, - sin 6*, a6l)}, 

and compute the optimal slopes with respect to them numerically. 

The first functional is the "average ropelength" which is the smallest length of a rope with unit 
thickness which is needed to make one twist of the double helix. It is a dual notion of thickness^ 
which is the radius of the "fattest tube" around a curve, or equivalently, the minimum global radius of 
curvature^. To be precise, the thickness of a curve F, which we denote by p(r), is given by the supremum 
of the radii r > so that a tubular r- neighbourhood of T (denoted by Nr{r)) can be embedded, i.e. the 
supremum of r > so that the normal discs to F of radius r at distinct points never intersect [3J IH [S] . 

The ropelength is the ratio of the length L{T) of F and the thickness p(F)[5J[3]. A knot that minimizes 
the ropelength in its isotopy class is called an ideal knot after Stasiak. In other words, an ideal knot is a 
knot whose fattest tubular neighbourhood maximizes a scale invariant functional on the space of knots 
defined by 

Vol(jVp(K)(i^)) 
lJ{K) 

in its isotopy class. The ideal double helix was computed by Sylwester Przybyl and Piotr Pierahski[7l[8]. 

The second is the "packing proportion" , in other words, the ratio of the volume of the fattest tubes 
and that of the circumscribed cylinder of the tubes (Figure [3|) . The best volume packing double helix 
was computed by Kasper Olsen and Jakob Bohr[9]. 

The last is an "average mutual energy^ , which is a central topic of the paper. An energy of knots 
is a functional on the space of knots which was introduced to produce a representative embedding for 
each isotopy class as an energy minimizer in the isotopy class. The first example was obtained as the 
renormalization of r~^-modified potential energy of knots pi3]. This energy is sometimes called the Mobius 
energy because it is invariant under Mobius transformations of the ambient space, but in fact, there are 
many other energies which are also invariant under Mobius transformations [TTl [T^. The energy blows 
up if a knot degenerates to a singular knot with double points. Therefore a knot can be deformed to 



'^to appear in Progress of Theoretical Physics Supplement, proceedings of "Statistical physics and topology of polymers 
with ramifications to structure and function of DNA and proteins" , Kyoto 2010. 
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decrease its energy without having self-crossing, and hence one can expect that the isotopy class would 
be kept the same during the deformation. 

Many kinds of energies have been studied intensively for about twenty years, forming a branch of 
mathematics called "physical knot theory" . In this paper we use two one parameter families of energies 
for ri(a)ur2(a), 

/ / dxdy {a > 1) and [ [ ~ exp{~kr) dxdy (fc > 0) , 

JTi{a)JT2{a) ^Fi (a) J Fa (a) 

where r denotes the distance \x — y\ between a pair of points. Namely, we consider the integration of 
modified Coulomb's potential and the screened Coulomb's potential (Yukawa potential). As both diverges 
because of non-compactness of the domain of integration, we take the "average per one twist" , in other 
words, we consider the integration on Ii{a) x T2{a), where Ii{a) is a subarc of ri(a) that winds up 
just once around the axis. As we only consider the interaction between two strands Ti{a) and r2(a), 
the energies are well-defined without renormalization. We compute the slopes that minimize the above 
energies and see how they depend on the parameters a and k. 

We will compare our numerical results with experimental data of DNA reported by Stasiak and 
Maddocks. 

We remark that our strategy of looking for optima double helices with respect to one of the three 
kinds of functionals above mentioned does not work for single helices. In fact, our functionals are all 
optimized when the slope a goes to -|-oo, i.e. when the curve approaches a straight line. 



2 An "ideal" double helix 

2.1 Thickness 

Let r be a curve. Let d : F x F R be the distance function: d{x,y) = |a; — y\. The doubly critical self 
distance of F[13l[3], denoted by dcsd(F), is the smallest positive critical value of d, i.e. 

dcsd(F) = min{|a; - y\ : t^, ty ± {y - x)}, 

where tx denotes the tangent vector to F at a;. Then the thickness of a knot F is given by [3] 

p(F) = min |min{radius of curvature at a;}, — dcsd(F) 1 . (2-1) 
[ X 2 J 

Remark that, for a fixed point xq, the radius of curvature at xq is a local quantity whereas miny^ajgHxo — 
y\ : txg,ty _L {y — xq)} is a global quantity. 

We remark that, in the case of a single helix with slope a, Maritan, Micheletti, Trovato, and 
Banavar[14) showed that there is a critical slope oo ~ 0.399805 ~ 2.512/27r such that the thickness 
is given by the half of the doubly critical self distance if a < oq and by the minimum radius of curvature 
if a > flQ. 



2.2 Average ropelength and the slope of an ideal double helix 

Let AL{a) be the length of F(a) per one twist. It is equal to An^/l + . Define the average ropelength 
of a double helix F(a) with slope a, denoted by ARL{a), by 

ARL{a) = — — = — , 

p{a) p[a) 

where p{a) denotes the thickness of r(a). Let us compute the slope a of an ideal double helix, i.e. the 
slope a that attains the minimum average ropelength. 

The thickness p{a) of F(a) is given as follows. The radius of curvature is constantly equal to 1 + a?. On 
the other hand, the doubly critical self distance is equal to the distance between a point (1, 0, 0) e Fi(a) 
and F2(a). The latter is given by 2 if a > 1 and by 

^J (cos Oq + 1)2 + sin^ Oq + g^Oq^ = \J 2 + 2coseQ + o^Oq^ , 
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where 6*0 {0 < 9o < n) is given by sin 6*0 = cl^Oq if a < 1. Therefore, by (|2.1I) . the thickness is given by 
the half of the doubly critical self distance for any a: 



p{a) 



1 (a > 1) 

W2 + 2 cos 6*0 + 026*0^ , (0 < 6lo < tt, sin^o = a^eio) (0 < a < 1). 



(2.2) 



When a > 1 we have ARL{a) — Att\/T+c? which takes the minimum value A\/2tt w 17.7715 at a = 1. 
The numerical computation using Maple implies that the average ropelength takes the minimum which 
is approximately equal to 17.003 at a « 0.82074 (Figure[T]), when the thickness p is 0.95614 and the ratio 
of the pitch and the thickness is P/p — 2TTa/p{a) ~ 5.3934. 




Figure 1: The average ropelength of a double 
helix with slope a 




Figure 2: Optimal tubular double helix with respect 
to the ropelength. This ideal double helix had al- 
ready been reported by Pieranski[8 , Fig. 12 (b). 



The above result is not new: This double helix with a straight axis that is optimal with respect to 
the ropelength (Figure [2]) had already been reported by Pieranski[8], Fig. 12 (b). 

The equation (j2.2p implies that the closest point in T2{a) from a point in ri(a) can be found in the 
opposite side when a > 1 and in an upward (or downward) direction when a < 1. It means that the curve 
of the contact of two fattest tubes around ri(a) and r2(a) is the straight axis when a > 1 i.e. when the 
ratio of pitch and thickness satisfies P/p> 2tt, and a helix between ri(a) and r2(a) when < a < 1 i.e. 
when P/p < 27r, as was explained by Stasiak and Maddocks|T5]. 



3 Best volume packing double helix 

Next we compute the slope of "best packed" tubular double helix, namely, a double helix that maximizes 
the ratio of the volume of the fattest tubular neighbourhood and the volume of the circumscribed cylinder 
of the tubes (Figure [3]) . This "packing proportion function" is given by 

TT p'^ia) ■ AL{a) _ 2VT 
" ^(l+p(a))2.27ra " a(l + p(a))2 ' 

where p{a) is the thickness of r(a). When a > 1 we have PR{a) = \/l + a? /2a which takes the maximum 
value at \/2/2 « 0.707107 at a — 1. The numerical computation using Maple implies that the packing 
proportion takes the maximum which is approximately equal to 0.7694 at a « 0.635805 (Figure|4]), when 
the thickness p is 0.832589 and the ratio of the pitch and the thickness is P/p ^ 4.7981. 

The best volume packing double helix had already been studied by Kasper Olsen and Jakob Bohr [5]. 
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Figure 3: Best volume packing tubular double helix and Figure 4: The packing proportion of a 
the circumscribed cylinder double helix with slope a 





Figure 6: The same one seen from an axial 
direction 



Figure 5: Best volume packing tubular double helix 
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4 Energy minimizing double helix 

4.1 r "-modified potential energy 
4.1.1 Mutual energy 

Define the r^"' -modified Coulomb's potential of ^2(0.) at a point a;o in ri(a) by 

dy 



r2(a) \xo-y\' 



de 



00 {{cos0 + lf + sm'^e + a^0^}' 



We assume a > 1 so that V^"\a) is finite. As a increases y(")(a) is dominated by the contribution 
of near-by points. Remark that it does not depend on the point xq € ri(a) because of the symmetry. 
Define the average cross-term r~°' -modified Coulomb's potential energy oi T (a) by 



A£;(")(a) = -AL{a) ■ V'^°'\a) = 27r(l + a^) 

2 ./-oo |(cosf 



d9 



sm 



The asymptotic behavior of AE'-"^ (a) as a goes to 00 or is given as follows: AS(")(a) - Coo ■ a as 
a +00 and AE^'^^ia) > Cq ■ -4=n: as a ^ +0, where Coo = 2'^-°'Tr (1 + t^)-! dt and Co = 
27r2""/_^(l+t^)-t<it. 

Numerical experiments using Maple imply that, for each a > 1, AE'-°'^{a) is a convex function of a. 
The graphs of AE^^^a) when a — 2 and 1000 are illustrated in Figures [7] and IH 
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Figure 7: (^-j 



Figure 8: AE^^^^^^a) 



Let a(a) denote the slope that minimizes AE^"'^. Their values obtained by numerical computations 
using Maple are shown in Table I. They are illustrated in dotted diamonds in Figure [9l The computations 
when a equals 1.5, 200, 400, and 1000 were first done by Miyuki Tani. One of the candidates for a fitting 
curve which is drawn by a blue curve in Figure ^ is given by 

14.0806074006923936 
■ ^ ia + 3.31656608512330786)2 ' 



4.1.2 Total energy when a = 2 

The energy AE^"^{a) is based on the interaction of two different strands Ti{a) and T2{a). If we further 
take into account the interaction between the same strand, we obtain the "total energy" . As the integral 
diverges at the diagonal set, we need renormalization to get a well-defined functional. 
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a 


a{a) 


P/p 


1.05 


1.01015 ^ 1.01018 


6.3471 


1.1 


l.Uzlzo ~ l.UzloU 


0.4159 


1.2 


l.U4oU5 ^ l.U4oll 


c c c: on 

6.5539 


1.5 


1.109545 ^ 1.109555 


6.9715 


1.75 


1 1/^OOn 1 "I/^OOO 


7.3124 


2 


1.21513 ~ 1.21515 


7.6349 


2.5 


1.30566 '-^ 1.30568 


8.2038 


3 


1.378780 ^ 1.378786 


8.6631 


4 


1.48104 ~ 1.4ollz 


n on c n 

9.3059 


5 


1.54403 ~ 1.54406 


9.7015 


D 


1.58435 ^ 1.58445 


9.9551 


o 

8 


1.63085 ~ 1.63095 


10.247 


1 n 
iU 


i.6oo96 ~ i.6oo98 


iU.4Uo 


15 


1.68555 - 1.68565 


10.591 


20 


1.69878 - 1.69881 


10.674 


40 


1.71655 - 1.71660 


10.786 


100 


1.72606 - 1.72617 


10.846 


150 


1.72810 - 1.72818 


10.858 


200 


1.72910 - 1.72913 


10.865 


400 


1.73056 - 1.73062 


10.874 


1000 


1.73139 - 1.73140 


10.878 



Table I. AE^"^^ -minimizing slopes a{a). The 
thickness is always 1 as a > 1. 




10 20 30 40 



Figm-e 9: The data of a{a) (diamond) and a fit- 
ting curve /(a) (blue) 



Figure 10: Fattest tubes around an energy minimizing double helix T{a). (left) a ~ 1 i.e. P/p — 27r. 
We conjecture that it is the limit of v4i?("-'-minimizer as a goes down to 1, and at the same time, the 
limit of i?ixp-minimizer as k goes down to 0. (center) a = 1.21514 i.e. P/p = 7.6349. Optimal with 
respect to AE'^°'^ when a = 2. (right) a = 1.7314 i.e. P/p = 10.878. We conjecture that it is the limit of 
Ai?("^-minimizer as a goes to -|-oo, and at the same time, the limit of i?ixp-minimizer as k goes to -l-oo. 
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Define the renormalized r ^ -modified Coulomb's potential at a point xq in Ti(a) bv|10) 



VtolM = 1™ ( / 



dy 



dz 



where Ng{xo) denotes the e-neighbourhood of Xq in ri(a), {y G ^i{a) : d{xo,y) < s}. Here d{xo,y) is 
the distance between ccq and y, either the arc-length or \xo — y\, both giving the same result|TB]. Define 
the average renormalized r^^ -modified Coulomb's potential energy of r(a) per one twist by 

Rob Kusner and John Sullivan showed that it takes the minimum value at a « 1.454 [S]. Remark that it 
is greater than the value of the optimal slope for AE^^\ This is because the "self energy" of Ti{a) is a 
decreasing function of a. 



4.2 ^ exp(— fcr)-energy 

Let Xq be a point in Ti(a). Define 

1 



r2(Q) 



, exp ( - fcV2 + 2 cos 61 + a^O^ ) , 

exp(-kr)dy^2 ^ , '-Vl + ^de, 

' ^ I V2 + 2cos(? + a202 



where r — |xo — y|. We assume /c > so that Vcxp{a) is finite. As k increases T4^xp(a) is dominated by 
the contribution of near-by points. Remark that it does not depend on the point xq e ri(a) because of 
the symmetry. Define 



/,\ L ii\ n ; exp ( — k^/2 + 2 cos 9 + ) 

Etiia) = -AL{a) ■ (a) = 47r(l + a') < ^ ^ ^ ' 



\/2 + 2cos0 + a^9^ 



de. 



Numerical experiments using Maple imply that, for each k > Q, Ei^-p{a) is a convex function of a. 
The graphs of Ei%{a) when fc = 0.1 and 200 are illustrated in Figures [TT] and [12] 



180- 



160 



140- 



120- 



100- 



80- 



60- 



1.4x 10" 



1.2x 10- 



1.x 10- 



6.x 10- "'L 



0.5 



1.5 



2.5 



1.5 



2.5 



Figure 11: Etp\a) 



Figure 12: Ei^''\a) 



Let acxp(fc) denote the slope that minimizes i?ixp. Their values obtained by numerical computations 
using Maple are shown in Table IL They are illustrated in dotted diamonds in Figure [T31 One of the 
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candidates for a fitting curve which is drawn by a red curve in Figure [13] is given by 

0.899179433232235592 

g(k) = 1.7314 3 . 

(fc + 1.17422998242145948)2 



k 




P/p 




0.01 


1.01947 


6.4055 




0.1 


1.10680 


6.9542 


1.7- 


0.3 


1.20560 


7.5750 




0.5 


1.29850 


8.1587 


1.6- 


0.75 


1.39250 


8.7493 




1 


1.45953 


9.1705 


1.5- 


1.5 


1.5409 


9.6818 


1.4- J 


2 


1.58628 


9.9669 


3 


1.6340 


10.267 


1.3-i. 


4 


1.65845 


10.420 


5 


1.6732 


10.513 


1.2-'> 


6 


1.6831 


10.575 


8 


1.69543 


10.653 


1.1-> 


10 


1.70282 


10.699 


■ > 


20 


1.7175 


10.791 


1.0- 


30 


1.72238 


10.822 





100 


1.72916 


10.865 




200 


1.7306 


10.874 





Figure 13: The data of aoxp(fc) (diamond) and a 
Table II. £;i^|,-niinimizing slopes aexp(fc>. The ^^^ting curve ff(/c) (red) 
thickness is always 1 as a > 1. 



4.3 Conjectures 

We like to end the section with conjectures implied by numerical experiments. 
Conjecture. Let a > 1 and fc > 0. 

(1) For fixed a, AE^'^^a) is a convex function of a. 
For fixed k, £'ixp(a) is a convex function of a. 

(2) The slope a{a) that minimizes ^i?'"^ (a) is an increasing function of a. 
The slope acxp(fc) that minimizes E'ixp is an increasing function of k. 

(3) Both {a(a) I a > 1} and {acxp(fc) I fc > 0} are bounded subsets of R, and we have lim a{a) — 

a—fl+O 

lim acxn(fc) = 1 and lim a(a) = lim acKuik) < \/3. 

5 Conclusion 

We have studied three kinds of functionals for double helices with straight axes and compute the optimal 
slopes with respect to them numerically (i.e. without mathematical proofs). 

If we consider the average ropelength, or equivalently, the ratio of the volume of the fattest tubular 
neighbourhood and the cube of the length, the maximum is attained by an ideal double helix, whose 
slope is 0.82074 (Figure [H c). 
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If we consider the ratio of the volume of the fattest tubular neighbourhood and that of a circumscribed 
cylinder of the tubes, the maximum is attained by a best volume packing double helix, whose slope is 
0.635805 (Figure [H b). 

If we consider two kinds of one parameter families of interaction energies between two strands of types 
r^" (modified Coulomb potential) and ^exp{—kr) (screened Coulomb's potential or Yukawa potential), 
the optimal slopes seem to be between 1 and \/3 in both cases, approaching 1 (or 1.73144) (Figure [Ml 
e,f) as the energy is dominated by the global (or respectively local) contribution, i.e. the contribution of 
near-by points gets weaker (or respectively stronger) . 

Stasiak and Maddocks[l5] reported that, in the case of DNA, the ratio P/p is 6.03 (Figure [141 d), 
which is within 4% error of that of optimal double helix if we use "r~-'^-potential energy", i.e. the limit 
of either interaction energy as the parameter approaches to the critical value when the energy fails to be 
well-defined, when the contribution of near-by points gets weakest. 




> a 



Figure 14: The slope a and the ratio of the pitch and the thickness P/p, and the optimal ones with 
respect to various functionals: (6) volume packing ratio (c) ropelength (ideal double helix) (d) DNA (e) 
the limit of the interaction energy AE^°''' as a \ 1 or Ecxp as fc \ (conjecture) (/) the limit of the 
interaction energy AE^"^ as a -l-oo or i?ixp as fc +oo (conjecture). Any point between e and / can 
be realized as an optimal slope with respect to some interaction energy. 
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